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A theory of charge transport in semiconductors in the presence of basal stacking faults is developed. It is 
shown that the presence of basal stacking faults leads to anisotropy in carrier transport. The theory is applied 
to carrier transport in non-polar GaN films consisting of a large number BSFs, and the result is compared 
with experimental data. 



8 The III-V nitride semiconductors and related com- 

9 pounds have attracted immense attention for optoelec- 

10 tronic devices 1 ' 2 covering a wide range of the electromag- 

11 netic spectrum as well as high-speed, high temperature 

12 electronic devices. For c-plane grown nitrides, built-in 

13 polarization field, although advantageous for two dimen- 

14 sional electron gas (2DEG) formation in transistors 3 , re- 

15 duces the overlap of electron-hole wavefunction, reducing 
is the oscillator strength for radiative transitions. More- 
17 over, c-plane based enhancement mode transistors suffer 
is from low threshold voltages 4 (Vth ~ 1 V) impeding their 

19 applications in safe circuit operation. A potential way 

20 to eliminate these effects is to grow nitrides in non-polar 

21 directions. However, heteroepitaxial, non polar and semi- 

22 polar nitrides films grown on lattice mismatched sub- 

23 strates contain «sf ~ 1-10 x 10 5 /cm basal stacking faults 

24 (BSFs) parallel to the [0001] direction 5-7 . Directionally 

25 dependent transport measurements in highly faulted non- 
26 polar nitride films show strong mobility anisotropy for 

27 both electron and hole transport, with mobility parallel 

28 to the [0001] direction significantly lower relative to the 

29 in-plane [1120] mobility 8-10 . This anisotropy has been 

30 associated with scattering from BSFs qualitatively, al- 

31 though no quantitative model exists. In this letter, we dc- 

32 velop a theory of carrier scattering phenomena from BSFs 

33 and a quantitative estimate of the transport anisotropy 

34 is presented with comparison to experimental data. 

36 

37 Stacking faults (SF) are 2D structural defects associ- 

38 ated with heteroepitaxial growth on lattice-mismatched 

39 substrates. Most prevalent SFs in wurtzite (WZ) nitrides 

40 is of the Ii-type requiring lowest energy of formation 11 . 

41 Scattering-contrast transmission microscope (TEM) 

42 imaging revealed 5 that the SFs in GaN are primar- 

43 ily of the Ti-type which corresponds to a stacking se- 

44 quence of (0001) basal plane • ■ ■ ABABABCBCBC 

45 A BSF can be thought as a thin zincblende (ZB) layer 

46 (up to three monolayers thick) perfectly inserted in the 

47 WZ matrix without broken bonds as shown in Fig. la). 

48 The built-in polarization difference between ZB and WZ 

49 structures will result in bound sheet charges ±ov at each 
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50 interface of the BSF layer. Consequently, the band edge 

51 will bend near the BSF. Fig. lb) shows a typical conduc- 

52 tion band diagram around a BSF for a n-doped GaN 

53 film along with the charges that are formed in the di- 

54 rection perpendicular to the BSF plane. The bending 

55 of conduction band edge inside the BSF is due the elec- 

56 trie field resulting from the polarization bound charge, 

57 whereas, bands bend outside the BSF due to the accu- 

58 mulation and depletion of mobile charges. We approxi- 

59 mate the accumulation charge as a sheet charge of den- 

60 sity n s at the centroid (t) of the charge distribution, i.e. 

61 p n {x) = en s 5(x — x ), where x = d + t, d being the 

62 width of a BSF and e is the electron charge. If Xd is 

63 the width of the depletion region, charge neutrality re- 

64 quires n 8 = XdN d , where is the donor density. Energy 

65 conservation across the BSF leads to the relation 

eN d x 2 d en s t _ /eoy en s\ d (1 \ 
2e s 2e s ~ V 2e s 2e s ) ' U 

66 where, e s is the permittivity of the semiconductor. In the 

67 limit Xd 3> (t, d), the depletion length is Xd — \J u^d/Nd- 
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FIG. 1. a) structure of basal stacking faults , b) charges across 
the basal stacking fault, c) schematic diagram of electric field, 
d) conduction band diagram including band-offset, e) delta- 
function model of barrier and QW. 
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The depletion width Xd increases with decreasing donor 
density and there exists a critical donor density N^ r for 
which, Xd becomes equal to the distance between two 
stacking faults. As a result the whole channel is depleted, 
and conduction ceases along the x direction. The critical 
donor density above which onset of conduction along x 
direction occurs is given by N^ r > a^dn^p. For a typical 
fault density iisf = 10 5 /cm and d = 0.8 nm for GaN 12 , 
the critical doping density is N^ r = 10 16 /cm 3 . In the rest 
of the paper, we will work in the regime where Nd N% r . 

In addition to the polarization charges at interfaces, 
the conduction band offset (AE C ) between ZB and WZ 
structures leads to a quantum well (QW) in the fault re- 
gion as shown in Fig. lb). For an applied bias, electrons 
tunnel through the barrier (in the depletion region) and 
QW (in the fault region) and then diffuse in the space 
between two consecutive stacking faults. To model the 
transmission coefficient of tunneling, we approximate the 
barrier and the QW as two delta functions (as shown in 
Fig.le) of strengths Si = eV^Xd and 5 2 = (AE C + eV v )d; 
where — a n d/2e s . The energy dependent coefficient of 
transmission through a single delta function potential is 
an analytically solvable problem 13 . If T tr ,i(e) and Tt rj2 (e) 
are the transmission coefficient from the barrier and the 
QW respectively, then the total coefficient of transmis- 
sion across the BSF is given by T tr (e) = T tT . t x(e)T tr! 2(e): 



conductivity across the BSF is then given by 



Ttr(e) = 



1 
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95 where, e is the energy of the incoming electron and m* is 

96 the effective mass at the band edge. In SF-free structures, 

97 the elemental current component along the x direction 

98 for an electron with velocity v k is given by j x = ev^Sfk', 

99 where Sf k is the modification of equilibrium Fermi-Dirac 

100 distribution caused by an applied electric field F app i. In 

101 the presence of SFs, a part (given by Eq. 2) of these 

102 carriers is transmitted through the barrier and the QW 

103 leading to an effective elemental current density jt = 

104 T tr {e)j k . Then, the total current density in the presence 

105 of SFs is given by 



J x = 2c 



d 3 k 
(2tt)3 



T tr (eH5f k , 



(3) 



we where, v k is the x component of electron group velocity 

107 in the state |r, k) and the factor 2 takes spin degener- 

108 acy into account. For a small applied electric field along 

109 the x direction, under the relaxation time approximation 
no (RTA), Sfk « {-df /de)r(k)v%F appl , where, f is the 
in equilibrium Fermi-Dirac distribution function and r{k) 
n2 is the momentum relaxation time due to scattering from 
in impurities, phonons etc present in the material. Inte- 
rn grating Eq.3 over all k-space and defining conductivity 
us by the relation J x — a(n,T)F ap pi, one obtains the car- 
n6 rier concentration (n) and temperature (T) dependent 
n7 conductivity in the presence of BSFs. The expression for 



a xx (T,n) 



3ne 2 J dee 1 / 2 cosh~ 2 (e/2k B T)T tr (e)T(e)(v^ 2 
~8tT /dee 3 / 2 cosh~ 2 (e/2fc B T) 

(4) 

where T is the equilibrium temperature, k B is Boltzmann 
constant. Since the BSF does not break the periodic 
symmetry in the y and z directions, inserting T tr (e) = 1 
for all energies, a similar expression can be obtained 
for conductivity (<J yv ) along the y direction. For hole 
transport, AE C should be replaced by valence band 
offset (AE V ) in the strength of delta function potential 
52, and a hole effective mass should be used instead of 
the electron effective mass. 

Equation (4) is the central finding of this work as 

30 it allows us to calculate experimentally measurable 

31 quantities such as conductivity and mobility in the 

32 presence of BSFs. As an application of the formalism 

33 constructed above we investigate charge transport in 
p-type m-plane (1100) GaN in the presence of BSFs, and 
compare the results with previous reported experimental 
data. Inclusions of different scattering mechanisms are 
necessary to evaluate the energy dependent momentum 
relaxation time r(e) appearing in Eq.4. Due to the 
high activation energy of acceptors ( Mg for p-type 
GaN has activation energy ~ 174 meV 8 ), a large 
fraction of the dopants remain neutral even at room 

42 temperature, acting as neutral impurity scatterers to 

43 hole transport. The fact that a high doping density 
is required to achieve appreciable hole concentrations 
results in high density of neutral impurities (NI) making 
neutral impurity scattering dominant even at room 

47 temperature. Accounting for the zero-order phase 

48 shift 14 , the momentum relaxation time of neutral 

49 impurity scattering is r^j = 20NHao/m*; where N is 

50 the density of neutral impurities and ao = 47re s 7i 2 /m*e 2 
is the effective Bohr radius. Momentum relaxation 

52 time due to ionized impurity scattering (Ti mp ) is 
calculated following the method of Brooks-Herring 
(BH) 15 . The fact that T imp ~ n~^ p ; where n imp is the 

55 density of ionized impurities, makes ionized impurity 

56 scattering important at high carrier concentrations 

57 (ionized dopant concentrations). For electron-optical 

58 phonon momentum- relaxation time (r p h), only phonon 

59 absorption has been considered due to the high optical 
phonon energy (E op = 0.092 eV > k B T) of GaN. The 
resultant momentum-relaxation time is calculated using 

62 Mathiessen's rule; i.e. r(e)^ 1 = t^) + r^ p + r~£. For 

63 numerical calculations, we assume the valence band 

64 offset A_E„=0.06 meV 11 . The hole effective mass of GaN 

65 is not well known and a wide range m* h = 0.4 — 2.4mo 
(toq is the rest mass of a bare electron) is found in 
existing literature 2 . In this work, we have assumed 

ea m\ ~ 1.8mo. Neutral impurity density is chosen to be 
N = 24 x p, where p is free hole concentration in GaN. 
This choice is not ad hoc; measurements 8 show that the 
neutral impurity density N = Na — P ~ 20 — 25 x p for 
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172 a wide range of acceptor densities Na- 

173 

174 Figure 2a) shows the calculated (solid lines) variation 

175 of a xx (along the c-axis) and a yy (parallel to a axis) as a 

176 function of hole density at room temperature juxtaposed 

177 with experimental 8 data (filled circles). It is apparent 

178 from the figure that the presence of BSF causes an 

179 appreciable reduction of a xx compared to a yy , resulting 
wo in anisotropic hole transport . The corresponding hole 

181 mobility pa = au/pe; i — x,y, is shown in Fig. 2b). The 

182 anisotropy in hole mobility increases with decreasing 

183 hole density. This stems out from the fact that, in the 

— 1/2 

184 limit p — > 0, the depletion width increases [xd ~ N A ) 

185 and one approaches the critical acceptor density limit 

186 N A r . As Xd increases, the strength of the barrier S\ 

187 increases, and T tr (e) — > 0. This results in vanishing /i xx 
las (a xx ; see Fig.2). 



wo We next investigate the electron transport in n-doped 

191 samples in the presence of BSFs. Experiment 8 shows 

192 that electron mobility in m-plane GaN decreases with de- 

193 creasing electron concentration. This phenomena is well 

194 described by charged dislocation scattering as pointed 

195 out by Weimann et. al 16 . The large lattice mismatch 
we between GaN and the substrate (for GaN on SiC, lat- 

197 tice mismatch ~ 2.9%) leads to high dislocation den- 

198 sities. For m-plane GaN grown on foreign substrates, 

199 TEM images reveal 17 ' 18 edge dislocation lines perpendic- 

200 ular to the m-plane ( along the [1100] direction) with 

201 non zero Burger's vector along [0001] direction. Each 

202 dislocation line acts as an acceptor-like trap with a line 

203 charge density pl — e/a, where a is the a-lattice con- 

204 stant of GaN. The electrons moving in the perpendicular 

205 plane to the dislocation line effectively feel a Coulomb 



206 potential V(r) = (pL/2ire s )K (r / 'Lp), where K (- ■ •) 

207 is the zeroth order Bessel function of 2nd kind and 

208 Ld = \JkBTe s /e 2 n' is the Debye screening length. Both 

209 free electrons, and bound electrons are taken account in 

210 the effective screening concentration (n 1 ) using Brook's 19 
2n formula n' — 2n — n 2 /Njj; where n is the free carrier con- 

212 centration. The momentum relaxation time (rdis) due 

213 to dislocation lines has been calculated by several au- 

214 thors for three 16,20,21 and two-dimensional electron gas 22 

215 in GaN. In the presence of BSFs, the x component of con- 

216 ductivity at room temperature due to dislocation scatter- 

217 ing can be calculated as 



Ve 2 / 47rVm*n dl sL D 



(5) 



219 where, Udis is the dislocation density and T(n. AE C ) is a 

220 dimensionless integral. The complex nature of the of the 

221 integral X(n, AE C ) prevents an analytical evaluation and 

222 one must rely on numerical techniques. For AE C = 0.27 

223 eV 11 , and in the non-degenerate limit, I(n,0.27) ~ 6, 

224 and /j, xx <~ \fn. Fig. 3a) shows the room temperature 

225 drift mobility with electron density for dislocation den- 

226 sity ridis — 2 x 10 9 /cm 2 . It not only captures the mobil- 

227 ity anisotropy due to BSF, but also the mobility varia- 

228 tion with carrier concentration; the behavior is similar to 

229 those measured experimentally 8,16 . For calculations, we 

230 used to*, 

231 relaxation time r" 1 - < dts -r , NI ^ , lmp ^ 

232 ered. Unlike p-doped GaN, the unintentionally doped C 

233 atoms serve as neutral impurities with a measured den- 

234 sity N ~ 5x 10 16 /cm 28 . Measurements of ref. 7 show free 

235 electron concentration using Hall technique is lower than 

236 actual dopant concentration in n-doped non-polar GaN 

237 samples. The difference between the dopant concentra- 

238 tion and the free electron concentration is attributed to 



m yy ~ 0-22too, and a resultant momentum 
rZl +tZI + rZl„ + is consid- 



4 



239 the acceptor traps associated with the dislocation lines. 

240 Moreover, the difference (Nd — n) ~ n ( u s varies from 

241 sample to sample indicating sample dependent density 

242 of dislocations. Assuming a single trap per a-lattice con- 

243 stant, the volume density of trapped charge is ndi S /a] 

244 and a charge balance equation can be written as 



a 1 + 2^ exp(E D /k B T) ' w 

245 where, Nc — 2(m*fcsT/27r?i 2 ) 3 / 2 is the effective density 

246 of states and Ed is the activation energy of the donor 

247 level. The screening effect on the donor activation 

248 energy is taken account by the empirical formula 

249 Ed = Edq — aNp 323 , where the screening parameter 

250 a = 2.1 X 10~ 5 meV.cm and E D o = 28 meV 23 . Using 

251 Eq. 6 and values of Nd and n from ref. 7 calculated dis- 

252 location densities are found to be in the range of n<jj S ~ 

253 1-10 x 10 8 /cm 2 . Electron mobilities is determined using 

254 the calculated values of sample-dependent dislocation 

255 densities and compared to experimentally measured 

256 values as shown in Fig. 3b). At low carrier densities, 

257 theoretical values are in agreement with experimental 

258 data; but they differ at high carrier concentrations. 

259 

260 An important observation is that the anisotropy 

261 in electron and hole mobilities decreases with in- 

262 creasing carrier concentration. At higher carrier 

263 densities, transmission of carriers occupying higher 

264 energy states approaches unity (T tr (e) — > 1) leading 

265 to /j, xx ~ Hyy. The corresponding carrier density is 

266 n (p ) ~ l/3ir 2 (m*AE c{v) d/h 2 ) 3 (~ 10 19 /cm 3 ), above 

267 which BSF-induced transport anisotropy vanishes in 

268 non-polar GaN. At high acceptor concentrations, the en- 

269 ergy levels of impurity atoms overlaps to form impurity 

270 bands. While conduction through this impurity band is 

271 important at low temperatures, at room temperature 

272 impurity band does not form an efficient transport 

273 path 2 and its exclusion from our model is quite 

274 justified. Also, the comparison between theoretically 

275 calculated drift mobility and experimentally measured 

276 Hall mobility is meaningful as at room temperature, the 

277 Hall factor rjj ~ l 23 . Our work presented here can be 

278 improved in three ways; i) at low carrier densities x c i is 

279 large; assumption of the delta function barrier breaks 

280 down and one should solve the Schroedinger equation 

281 numerically to obtain exact transmission coefficient 

282 T trl (e), and ii) at low doping levels (when is large) or 

283 for high stacking fault density, carriers are confined in a 

284 quasi-two dimensional space rather than moving in three 

285 dimensions. Hence, one should use a two-dimensional 

286 analogue of Eq.4 for an exact evaluation of \i yy , and iii) 

287 by incorporating the contributions of quasi-localized 12 

288 electrons in the shallow state of the thin QW located in 

289 the fault region. 

290 

291 High-quality (low defect densities) non-polar GaN is 



292 currently under development. It is apparent that with 

293 decreasing SF density, transport anisotropy will decrease 

294 and eventually vanishes for BSF free non-polar GaN 

295 films. The theory presented here is not only applicable 

296 to GaN, but to any semiconductor (including other 

297 members of the nitride family such as InN 6 and A1N 7 ) 

298 having BSFs with proper choice of material parameters. 

299 For example, the semiconductor zinc-oxide (ZnO) shows 

300 SFs similar to GaN with a„ — .057 C/m 2 and the band 

301 offsets 25 are AE C = 0.147 eV and AE V = 0.037 eV. 

302 

303 In summary, we have presented a transport theory in 

304 semiconductors containing a large number of basal stack- 

305 ing faults. The theory is applied to understand the exper- 

306 imentally observed transport anisotropy in m-plane GaN 

307 and a reasonable agreement between theory and experi- 

308 ment is obtained. Two critical limits of carrier concen- 

309 tration are derived, one where transport across a basal 

310 stacking fault ceases and another, where the presence of 
3n basal stacking faults can be ignored in the context of 

312 charge transport. 
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